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New  Data  Evaluation  Technique  for  Electron  Tunneling  Spectroscopy. 


Vladimii'  A.  Uki'aintsev 


University  of  Pittsbui’gh 
Department  of  Chemistry 
Surface  Science  Center 
Pittsbur’gh,  PA  15260 


A  systematic  study  of  Local  Density  Of  States  (LDOS)  deconvolution  from 
tip-surface  tunneling  specna  is  reported.  The  one-dimensional  WKB  approximation 
is  used  to  simulate  the  process.  A  new  technique  for  DOS  deconvolution  from  the 
electron  tunneling  spectroscopy  data  is  proposed.  The  differential  conductivity 
normalized  to  its  fit  to  the  turmeling  probability  function  is  used  as  a  method  of 
recovering  sample  DOS.  This  explicit  procedm'e  does  not  use  unconstrained 
parameters  and  reveals  a  better  DOS  deconvolution  in  comparison  with  other 
techniques.  The  advantage  of  tliis  new  method  is  its  feasibility  for  extracting  two 
important  physical  parameters  from  experimental  turmeling  spectr’a;  (i)  local  surface 
potential,  and  (ii)  tip-sample  distance.  These  values  are  the  par'ameters  used  in  the 
proposed  fitting  procedru'e.  The  local  surface  potential  and  the  tip-sample  distance 
retrieval  ar'e  demonstr'ated  by  means  of  numerical  simulations.  Compar'arive 
Scanning  Turmeling  Spectroscopy  (CSTS)  is  proposed  as  a  new  approach  to 
eliminate  the  influence  of  the  tip  condition  on  the  surface  LDOS  recovery. 


I.  Introduction 


The  first  obsei'vation  of  the  special  density  of  electi'onic  states  in  metal- 
oxide-superconductor  tunneling  experiments  was  made  by  Giaver.  ^  This  work  led 
to  a  reexamination  of  the  one-elecfron  theoiy  of  tunneling  that  has  dominated  the 
field  for  more  than  30  years^ .  The  new  many-body  ttansfer-Hamiltonian  approach 
in  tunneling  cuixent  calculation  was  pioneered  by  Bardeen^  and  developed  by 
Cohen  et  ,  Duke  et  a!.' ,  Appelbaum  and  Brinkman^ ,  Carol!  et  al.^  and  many 
others.*  This  foimalism  is  used  today  in  most  applications  of  the  tunneling  theoiy. 

The  invention  of  Scanning  Tunneling  Microscopy  (STM)  by  Binnig, 

Rohrer,  Gerber  and  Weibel^  stimulated  an  additional  interest  in  the  theoiy  of 
tunneling.  Using  the  transfer-Hamiltonian  foimalism^,  thi'ee-dimensional 
tunneling  in  the  STM  was  studied  by  Tersoff  and  Hamann^° ,  Garcia  et  al}\ 
Feuchtwang  et  andLang.^"* 

With  the  advent  of  STM,  a  Scanning  Tunneling  Spectroscopy  (STS)  with 
atomic  spatial  resolution  became  possible.  The  theoiy  of  STS  was  considered  by 
Feuchwang  et  al}^,  Selloni  et  and  Lang.^^  It  was  concluded  that  in  the 
general  case  of  three-dimensional  tunneling,  cmxent  cannot  be  calculated  as  a 
simple  convolution  of  sample  Densities  Of  States  (DOS)  and  tip  DOS  with  the 
“effective  matidx  element  for  tunneling”.  The  simple  relation  between  sample 
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Stroscio,  Feenstra  and  Fein*^  proposed  an  effective  solution  for  the  later 
problem.  It  was  found  that  by  normalization  of  the  differential  conductivity  by  the 
total  conductivity  one  can  effectively  remove  the  dependence  of  the  tunneling 
cmi'ent  on  the  tip-sample  distance.  The  result  was  qualitatively  generalized  and  it 
was  concluded  that  this  noimalization  procedure  will  cancel  out  both  tip-sample 
distance  and  tip-sample  voltage  dependencies  of  the  tuimeling  transmission 
probability  in  the  differential  conductivit\  .* '  It  was  claimed  that  the  (dL''dV)/(I/V) 
function  is  “a  relatively  direct  measure  of  the  surface  DOS”.^^  This  conclusion 
was  confirmed  by  Lang^^  in  thi'ee-dimensional  transfer-Hamiltonian  numerical 
calculations  for  tunneling  between  Na  and  Ca  single  atoms  positioned  on  infinite 
planar  metal  electi'odes.  The  noimalization  teclmique  was  also  justified  by 
Hamers. Simmon’s  foimulas^^  (WKB  approximation)  and  constant  DOS  for  the 
tip  were  used  to  simulate  the  tunneling  cuiient  dependence  on  the  tip-sample 
voltage.  The  recovered  (dI/dV)/(I/V)  function  has  shown  qualitative  agi’eement 
with  the  input  sample  DOS.  It  is  worthwhile  to  note  that  in  this  case  less 
comespondence  can  be  found  between  input  and  output  sample  densities  of  states 
in  comparison  with  the  “exact”  simulation  done  by  Lang.’^  It  is  also  important 
that  in  both  calculations^ metal-like  sample  DOS  were  used. 

Su'oscio  et  al}^  experimentally  tested  the  influence  of  the  tungsten  tip  DOS 
on  the  deconvoluted  sample  density  of  states.  The  DOS  of  a  Ni  sample, 
deconvoluted  by  the  (dI/dV)/(I/V)  technique,  shows  no  tip-related  features.  This 


conductivity  (for  details  see  Section  III.l).  Finally,  Feenstra  was  able  to  obtain 
reproducible  and  reliable  electronic  spectta  of  III-V  semiconductors.^^  However, 
with  all  these  empirical  adjustments,  it  is  increasingly  difficult  to  keep  explicit  the 
fundamental  meaning  of  the  noiinalization  procedure  as  well  as  to  justify  its 
ultimate  accuracy. 

In  this  paper  a  new  technique  for  DOS  deconvolution  from  elecU'on 
tunneling  specti'oscopy  data  is  introduced.  To  recover  the  sample  and/or  tip  DOS 
one  may  use  the  differential  conductivity'  normalized  by  its  jit  to  the  tunneling 
probability  function.  The  one-dimensional  and  semi-classical  WKB 
approximation  was  used  to  model  electi'on  tunneling  between  the  STM  tip  and  the 
sample.  As  shown  by  simulations,  in  all  studied  cases  this  technique  reveals  better 
deconvolution  of  the  input  DOS  in  comparison  with  any  existing  technique."  This 
procedui'e  does  not  use  unconstrained  parameters  except  the  analytical  foini  of  the 
tunneling  probability  function.  Since  one  should  employ  the  same  tunneling 
model  to  justify  any  other  DOS  deconvolution  technique,  the  suggested  procedure 
is  more  consti'ained  and  certainly  more  expUcit  than  others. 

The  advantage  of  this  new  method  is  the  feasibilitv'  of  exti'acting  rtvo 
important  physical  pai’ameters  from  experimental  tunneling  spectra;  (i)  “local  work 
function,”  and  (ii)  tip-sample  distance.  These  parameters  ai’e  used  in  the  proposed 
fitting  of  the  differential  conductivity  by  the  tumieling  probability  function.  The 
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II.l  Methods  of  Tunneling  Current  Calculation. 

The  tunneling  current  between  two  weakly  bounded  electi'odes  using  first- 
order  perturbation  theory  is  (refs.  6,  10,  12); 

IIS,  V)  =  \m,^,^‘s(e,  -  £j[f(£,  -  ev)  -  r(£j],  (1) 

^  t,S 

where  V  is  the  sample  bias  with  respect  to  the  tip;  5  is  the  tip-surface  separation, 
is  the  turmeling  matrix  element  between  states  if/j.  of  the  tip  and  y/,  of  the 

surface;  f{E)  is  the  Fermi-Duac  distribution  function;  and  E^  are  the  energies 
of  states  and  y/^,  respectively,  in  the  absence  of  turmeling.  Both  energies  are 
referenced  to  the  surface  Fermi  level.  Positive  curuent  indicates  electr'on  mnneling 
from  tip  to  surface. 

The  turmeling  matrix  element  can  be  calculated  as  follows  using  the  method 
of  Bardeen^: 


where  the  integral  is  calculated  over  an  arhitrary  surface  lying  entirely  within  the 
vacuum  region  separating  the  tip  and  the  surface. 
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The  tunneling  probability  is  a  strong  function  of  the  parallel  component  of 
the  energy  £||.  At  each  particular  total  energy  E,  the  DOS  with  zero  £|| 

component  is  heavily  weighted  by  the  tunneling  probability  in  equation  (3). 
Therefore,  p^is)  and  p^-is)  may  be  approximated  by  the  surface  and  the  tip 
densities  of  electronic  states  with  *  0  (Ref  18).  This  approximation  may  in 
pait  be  a  justification  for  ignoring  the  tunneling  electi'on  momentum  in  equation 

(3). 

One  can  simplify  equation  (3)  at  low  surface  temperature,  kT  «  eV  (k 
is  Boltzmann’s  constant),  by  using  the  step  function  instead  of  the  Feimi-Dirac 
distiibution  function.  The  density  of  the  tunneling  cuirent  is  tlien; 


2Ke  f  YY  /  \  f  \  f  \ 

J{S,  7)  =  —  —  I  t{s,  V,  e)p^{e)p.SE  -  eV)dt 


Assuming  E\i  »  G  for  the  electronic  states  which  influence  the  tunneling 
current,  the  fransmission  probability  can  be  written  in  the  following  form; 


Expressions  similar  to  (5,  6)  were  employed  in  most  cases  in  an  attempt  to 
describe  semi-quantitatively  the  tunneling  specti'oscopy  with  spatial 
resolution.*^’'^’^^’^^’^^  These  equations  were  obtained  for  the  weakly  bounded 
electrodes  in  the  one-dimensional  and  semi-classical  WKB  approximation  in  the 
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with  respect  to  the  tip  and  the  surface  densities  of  states.  One  can  accomplish  this 
by  substituting  ^  =  {e  -  ®%)  in  equations  5  and  6.  Then  the  differential 
conductivity  is; 


dl{S,  V)  ^  ^ 
dV  ^2 


+t'(s,  ^)p,(c  +  + 


C--V2 


eV/ 
~  /2 


dl 


1  t'(s,  c)p.(#  +  ^%) 


e  V  / 
”  /2 


a(#  -  - 

dp,U  -  ,J 

d# 


(8) 


where 


2  W. 

t’{s,  i)  ^  t(s,  V,  e)  =  expj  -2sJ-^  (O  - 


(9) 


is  the  appropriate  presentation  of  the  tunneling  transmission  probability  (6). 

Equation  (8)  is  symmetric  with  respect  to  the  surface  and  tip  DOS.  The 
r’(s,  is  a  probability  of  tunneling  to  or  from  the  tip  Fermi  level  and  the 
®%)  is  a  probability  of  tunneling  to  or  from  the  sample  Fermi  level.  As 
shown  in  Figure  1,  at  positive  V  (sample  bias  with  respect  to  tip) 
r'(s,  •%)  »  t’(s,-  ^%).  In  this  case  the  differential  conductivity  can  be 
crudely  estimated  as: 
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PsiO)  =  ^Ps(0)  =  0,  and  the  influence  of  the  tip  DOS  is  reduced  at  negative 
sample  biases  (11).  Therefore,  the  sample  density  of  occupied  electronic  states 
has  a  higher  chance  to  reveal  itself  at  negative  tip-sample  voltages  despite  the 
tunneling  probability  vanishing  (cf.  eqns.  8  and  1 1).  This  effect  will  be 
demonsti-ated  by  numerical  simulations  in  Section  IV. 

III.  Surface  and  Tip  DOS  Recovery  by  Normalization  to  Tunneling 
Probability. 

The  simplified  equations  (10,  11)  show  clearly  the  method  of  DOS 
deconvolution  by  using  the  differential  conductivity.  The  best  one  can  do  is  to 
honnalize  the  dl/dV  to  T  (s,  at  positive  sample  bias  and  to  t'{s-  ® %)  at 
negative  sample  bias.  Then  the  noimalized  differential  conductivity  will  be 
propoitional  to  the  sample  (10)  or  to  the  tip  (11)  DOS.  As  also  seen  from 
equations  (10,  11),  one  may  expect  an  influence  of  the  first  derivatives  of  the 
coiTesponding  DOS  on  the  recovered  specti’a.  Nevertheless,  all  existing  methods 
of  DOS  recoveiy  from  the  tuimeling  specttoscopy  data  ultimately  use  this 
approach.^^’^^’^’  The  difference  between  vaiious  techniques  is  in  the  different 
methods  of  recovering  the  tuimeling  transmission  probability. 

III.l  Methods  Based  on  Total  Conductivity  Utilization. 

.A.S  was  found  empirically  in  many  practical  cases'''^'  '^^’^  the  total 
tunneling  conductivity',  I/V,  follows  reasonably  well  the  desired  noraialization 
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probability  or  the  coiTesponding  differential  conductivity,  cTconsti^i  ^)-  Hence, 
normalization  of  the  differential  conductivity  to  the  total  conductivity  will  not 
recover  the  constant  input  DOS,  Ps{<^)  =  /9(.(^)  =  1. 


Significant  deviation  fi'om  the  v)  was  also  obtained  for  the 

modified  total  conducthity  (Figure  2A).  In  this  case  double  exponential  weighting 
of  the  broadened  and  extended  total  conductivity  was  used  to  approximate  the 
tunneling  transmission  probability  (Refs.  17,  23); 


X  ex 


p(-akl)dc, 


(14) 


where  a  is  the  best  parameter  obtained  duiing  the  |l(s,  \/)|  fitting  by  the 
exponential  function  exp(=!v|)  and  A  is  the  broadening  width.  In  all  simulations  A 
=  1  eV  was  used. 


Tunneling  cuiTent  versus  voltage  is  always  known  in  a  limited  voltage 
range.  Hence,  to  avoid  decay  of  the  convoluted  total  conductivity  in  the  proximity 
of  the  limits  one  should  somehow  extend  the  values  of  the  tunneling  current  or 
total  conductivity  far  beyond  the  limits  of  the  voltage  range.  For  variable  tip- 
sample  separation  tunneling  spectroscopy,  Feensfra^^  has  suggested  extending  the 
cun'ent  values  with  constant  value  equal  to  the  value  of  cument  at  the  largest 
measured  voltage  (the  extension  should  be  done  separately  for  positive  and 
negative  voltages).  For  constant  tip-sample  separation  spectroscopy  it  would  be 
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where  Bq  is  the  coefficient  used  at  constant  tip-sample  separation. 

DOS  deconvolution  accomplished  by  nonnalization  of  the  differential 
conducti'V’ity  to  the  modified  total  conductivity  (Ref.  23)  shows  a  systematic  error 
in  the  case  of  valuable  tip-sample  separation  (Figure  2B).  Therefore,  in  the 
framework  of  the  one-dimensional  WKB  simulations,  all  techniques  which  employ 
the  total  conductivity  for  DOS  recoveiy  reveal  significant  deviation  from  the  input 
tip  and  sample  constant  densities  of  states. 

IIL2  Method  Based  on  Differential  Conductivity  Fitting. 

The  differential  conductivity  at  a  constant  tip  and  surface  DOS  (eqn.  13) 
follows  closely  the  appropriate  tunneling  probability  at  positive  as  well  as  at 
negative  sample  biases  (Figure  1).  Therefore,  this  function  is  a  suitable  one  for 
DOS  recoveiy  by  normalization  (eqns.  10,  11). 

If  the  sample  Fermi  level  is  located  inside  the  band  gap  then 
p^iO)  «  PtiO)  and  the  tunneling  differential  conductivity  should  vanish  at 
negative  sample  bias  (eqn.  11).  The  symmetiic  tunneling  probability  function  (13) 
is  not  appropriate  for  the  fitting  of  this  asymmetiic  differential  conductivity.  The 
situation  can  be  fixed  by  employing  an  asymmetric  timneling  probability  function: 

f(s,  v)  =  Aj.t'(s,  +  AsT'{s-  ®%),  (17) 
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predicted  that  for  these  “up”  and  “down”  atoms  the  difference  in  nonnal 
coordinate  is  about  0.5  A.  Therefore,  for  the  STM  probe  placed  above  the  dimer 
the  density  of  states  of  the  “up”  atom  of  the  asymmetiic  dimer  are  enhanced  in  the 
tunneling  spectimm  by  a  factor  of  3  with  respect  to  the  density  of  states  of  the 
“down”  one.  In  this  analysis  the  decay  factor  has  been  assumed  to  be  k  -  1  A'^ 
(Ref  30 ). 

Therefore,  the  fitting  of  experimental  differential  conductivity  a  (f/  v)  by 

^consti^'  the  following  DOS  recoveiy  cannot  be  considered  to  be  a 

quantitative  procedure  but  rather  should  be  viewed  as  a  qualitative  one.  However, 
the  proposed  technique  re\  eals  the  best  DOS  recoveiy  at  least  in  the  frame  of  the 
WKB  simulations. 

IV.  Numerical  Simulation  of  the  Electron  Tunneling  Spectroscopy. 

WKB  approximation.  One-Dimensional  Case. 

In  the  simulations  presented  in  Figures  1-6,  the  following  values  were  used: 

B  =  12  eV**^^  coiTesponding  to  the  tip-surface  separation  of  -  11.7  A,  and  the 

average  tip-sample  work  function  0=3  eV.  Such  values  were  obtained  fi’om  the 

evaluation  of  tunneling  spectra  of  Si(001)-(2xl)  (Ref  3 1 ),  and  are  used  here  to 

better  model  the  experimental  simation. 

In  Figure  3  simulation  for  a  sine  tip  and  a  sine  surface  DOS  is  presented. 
The  periods  of  the  tip  and  the  sample  DOS  are  1  eV  and  0.5  eV,  respectri  ely.  As 


To  diminish  the  influence  of  tip  DOS  a  new  technique  teiTned  Comparative 
Scanning  Tunneling  Spectroscopy  (CSTS)  is  proposed.  The  effectiveness  of  the 
approach  is  demonstrated  in  Figure  6.  The  deconvoluted  DOS  shown  in  Figure  5A 
was  used  as  a  reference  spectinim.  Two  features  were  added  to  the  input  DOS  of 
the  reference  sample  on  both  sides  of  the  band  gap  in  order  to  get  a  probe  sample 
input  DOS  (Figure  5B).  The  same  tip  DOS  was  used  in  both  simulations.  The 
difference  of  input  sample  densities  of  states  (Fig.  6,  dashed  line)  as  well  as  tip 
DOS  and  output  comparative  spectra  are  presented  in  Figure  6. 

The  DOS  averaged  over  macroscopic  surface  area  can  be  obtained  by 
traditional  and  inverse  photoemission  spectroscopies.^^  Therefore,  as  soon  as  the 
reference  DOS  of  the  majority  surface  sites  is  known,  the  comparative  analysis 
provides  absolute  spectra  of  surface  elecuonic  states  with  high  spatial  resolution. 
To  get  all  advantages  of  the  CSTS,  both  tunneling  spectin  (reference  and  probe) 
should  be  obtained  with  exactly  the  same  tip.  STM  images  are  veiy  sensitive  to 
the  tip  condition  and,  hence,  the  later  can  be  precisely  monitored. 

Employment  of  the  tunneling  probability  function  (13)  for  the  comparative 
site  specific  tunneling  specmoscopy  on  Si(001)-(2xl)  revealed  excellent 
reproducibility  on  a  day-to-day  and  a  tip-to-tip  basis.^^  Application  of  this 
comparative  approach  in  the  case  of  dl/dV  nonnalization  by  the  modified  total 
conductivity  is  questionable,  partially  because  of  possible  changes  in  local  surface 
potential,  tip-surface  distance  and  effective  tip  area  from  one  surface  fearare  to 


time,  slow  convergence  and  larger  discrepancy  caused  by  high  coirelation  of  the 
fitting  parameters  (eqns.  9,  17)  was  observed  in  some  cases  of  highly  modulated 
tip  and  sample  densities  of  states.  Such  situations  never  occuired  during 
experimental  data  evaluation.^' 

However,  the  data  presented  in  Figure  7  can  be  considered  as  a  qualitative 
example  only.  One  should  clearly  realize  that  all  simulations  reported  here  were 
done  in  the  one-dimensional  and  semi-classical  WKB  approximation.  The  proper 
test  of  the  method  has  to  be  done  in  experiments  or  with  a  much  more 
comprehensive  theoretical  approach  (see  refs.  15,  33  ). 

V.  Conclusion 

A  systematic  study  of  surface  DOS  deconvolution  fi’om  the  tip-sample 
tunneling  electron  spectroscopy  data  was  undertaken.  The  one-dimensional,  semi- 
classical  WKB  approximation  was  used  to  simulate  the  electi'on  tunneling  between 
the  metal  tip  and  the  surface.  Using  the  simplest  case  of  constant  tip  and  sample 
densities  of  states  (Fig.  2),  it  is  shown  that  all  cuiTent  methods  of  DOS  recovery 
based  on  the  utilization  of  the  total  conductivity  have  a  systematic  eiTor.  In  this 
paper  a  new  technique  for  DOS  deconvolution  from  tunneling  spectroscopy  data  is 
introduced.  The  differential  conductivity  normalized  to  its  fit  to  the  tunneling 
probability  function  is  employed  as  a  method  of  recovering  sample  DOS.  As 
shown  by  numerical  simulations,  this  technique  reveals  a  better  deconvolution  of 
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possible  changes  in  local  surface  potential  from  one  surface  feature  to  another. 
The  influence  of  this  parameter  on  the  modified  total  conductivity  is  unknown. 
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Figure  3.  Tip  influence  on  DOS  deconvolution.  Top  graph:  input  sample  (upper) 
and  tip  (lower)  DOS.  Bottom  graph:  solid  line  is  the  DOS  deconvoluted  by  F 
normalization,  dotted  line  is  the  DOS  deconvoluted  by  noimalization  to  the 

modified  total  conducti\ity  at  constant  tip-sample  separation  ,  dash-dot  line  is 
the  DOS  deconvoluted  by  normalization  to  the  modified  total  conductivity  at 


variable  tip-sample  separation 


Values  used  in  the  simulation  are  the  same 


as  for  Figure  1. 


Figure  4.  Problem  of  sample  DOS  recovery  below  the  Feimi  level.  Top  graph: 
input  sample  (upper)  and  tip  (lower)  DOS.  Bottom  graph:  solid  line  is  the  DOS 
deconvoluted  by  F  normalization,  dotted  line  is  the  DOS  deconvoluted  by 
normalization  to  the  modified  total  conductivity  at  constant  tip-sample  separation 

^ ,  dash-dot  line  is  the  DOS  deconvoluted  by  noimalization  to  the  modified  total 

conductivity  at  vaiiable  tip-sample  separation  ^y(/.  Values  used  in  the  simulation 
are  the  same  as  for  Figure  1. 


Figure  5A.  Deconvolution  of  tire  simple  semiconductor  sample  DOS.  Top  graph: 
input  sample  (upper)  and  tip  (lower)  DOS.  Bottom  graph:  solid  line  is  the  DOS 
deconvoluted  by  F  noimalization,  dotted  line  is  the  DOS  deconvoluted  by 
normalization  to  the  modified  total  conductivity  at  constant  tip-sample  separation 
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variable  tip-sample  separation 
for  Figure  1. 


Values  used  in  the  simulation  are  the  same  as 


Figure  7A.  Work  function  (WF)  and  tip-sample  distance  (TSD)  recovery  in  the 
case  of  the  simple  semiconductor  sample  DOS.  The  tunneling  cuiTent  vs.  tip- 
sample  bias  dependencies  were  fitted  by  the  tunneling  probability  function  (17). 
Top  graph:  solid  line  is  the  dependence  of  recovered  tip-sample  distance  on  input 
work  function,  dotted  line  is  the  input  tip-sample  distance.  Bottom  graph:  solid 
line  is  the  dependence  of  the  recovered  work  function  on  input  work  function, 
dotted  line  is  the  input  work  function. 


Figure  7B.  Work  function  (WF)  and  tip-sample  distance  (TSD)  recoverv'  in  the 
case  of  the  simple  semiconductor  sample  DOS.  The  tunneling  cuixent  vs.  tip- 
sample  bias  dependencies  were  fitted  by  the  tunneling  probability  function  (17). 
Top  graph:  solid  line  is  the  dependence  of  the  recovered  work  function  on  input 
tip-sample  distance,  dotted  line  is  the  input  work  function.  Bottom  graph:  solid 
line  is  the  dependence  of  the  recovei'ed  work  function  on  input  work  function, 
dotted  line  is  the  input  work  function. 
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